Signatures of the classical transition state in atomic quantum spectra 
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We perform quantum mechanically exact calculations of resonances in the spectrum of the hydro- 
gen atom in crossed external fields and establish a close connection between the classical transition 
state in phase space and features in the quantum spectrum. By varying the external field strengths, 
structures are revealed which are surprisingly similar to the quantized energy levels of the classical 
electron motion in the vicinity of the saddle point obtained with an approximation of the potential. 
The results give clear evidence for signatures of the transition state in quantum spectra. 

PACS numbers: 32.60.-K, 32.80.Fb, 82.20.Db 



I. INTRODUCTION 



Transition state theory, which has its original appli- 
cations in the theory of chemical reactions (see, e.g., 
P, 0, can also be applied to many dynamical systems 
which evolve from an initial ( "reactants" ) to a final state 
("products"). The key concept is based on trajectories 
in the classical phase space of a dynamical system. In 
particular, regions of the phase space which identify the 
reactant side and those which represent the product side 
are introduced. Of special interest are trajectories that 
connect the reactant side with the product side, because 
they describe the "reaction." The theory postulates the 
existence of a minimal set of states which is passed by all 
of these "reactive" trajectories, and this set is called the 
transition state. The transition state has no intersection 
with a "nonreactive" trajectory, and therefore describes 
a boundary between reactants and products in the phase 
space. 

Among the large number of applications of the tran- 
sition state theory, it has also proven to be important 
for understanding the ionization mechanism of atoms. In 
the case of the hydrogen atom in an external electric field 
the "reaction" corresponds to the ionization. The poten- 
tial of the system is plotted in Fig. Q] (solid red grid) , in 
which the Stark saddle point is clearly visible. Classical 
trajectories describing an ionization must pass the saddle 
point, and therefore the transition state must be located 
in its vicinity. An early work by Clark et al. [J] can 
be regarded as a first step in discussing the ionization 
mechanism of the hydrogen atom in crossed electric and 
magnetic fields in this context. They used a quadratic 
expansion of the potential around the Stark saddle point 
(see Fig. [TJ dashed blue grid) to calculate quasi-bound 
states confined to the vicinity of the saddle. Further 
progress was achieved by Jaffe et al. [5, 6] who discussed 
transition state theory for systems without time-reversal 
symmetry, and found the transition state of the planar 
(i.e., two-dimensional) hydrogen atom in crossed exter- 
nal electric and magnetic fields. With the help of meth- 
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FIG. 1: (Color online) Potential of the hydrogen atom in 
an external electric field in the (£,t/)-plane (cf. Eq. ((2|) with 
£ = x — x s denoted by the solid red grid. Reactive trajectories 
must pass the Stark saddle point at (£ = 0, y — 0) in whose 
vicinity the transition state is located. The dashed blue grid 
describes a quadratic approximation around the saddle point 
(cf. Eq. ©). 



ods from nonlinear dynamics Jaffe et al. H Q discussed 
the chaotic ionization mechanism, and were able to ex- 
plain the existence of electrons ionizing promptly and 
electrons ionizing with delay after repeated encounters 
with the core region. The problem that the practical 
applicability of transition state theory was restricted to 
low-dimensional systems was overcome by an algorithmic 
procedure which allows identifying the transition state in 
higher-dimensional systems @, [a @, and is based on a 
normal-form representation of a power series expansion 
of the Hamiltonian. The quadratic approximation used 
by Clark et al. |4[ coincides with the lowest order in the 
normal form expansion of the Hamiltonian introduced by 
Uzer et al. Q- 

Even though the transition state for the hydrogen 
atom in crossed fields had been identified, the question 
of whether or not the classical trajectories in its vicin- 
ity leave a signature in the exact quantum spectrum re- 
mained unanswered since no such quantum spectra were 
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available. We have succeeded in calculating these spec- 
tra, and it is the purpose of this paper to establish the 
close connection between the classical motion of an elec- 
tron confined near the transition state and exact quan- 
tum resonances. To calculate the electron trajectories 
near the transition state we adhere to the simplest ap- 
proximation and use the power series expansion of the 
potential around the Stark saddle point up to second or- 
der introduced by Clark et al. Higher orders would 
rapidly lead to a drastic increase of the computational ef- 
fort, and, as will be seen from the results, the agreement 
of the quantized energy levels of the electron's motion 
with resonances in the quantum mechanically exact spec- 
trum is already extremely good for the simple expansion. 
Thus, the expansion is completely sufficient to reveal the 
signatures of the transition state in the quantum spec- 
trum. 

The crossed-fields hydrogen atom is an ideal candidate 
to search for the connection between the classical transi- 
tion state and its counterpart in a quantum system be- 
cause the lowest-order approximations to the transition 
state as well as the quantum states are accessible with 
reasonable numerical effort. Chemical reactions with 
their large number of degrees of freedom are far more 
complicated and do not allow for the insight we gain in 
the case of the hydrogen atom, where only three degrees 
of freedom exist. 

The paper is organized as follows. The system and the 
approximation around the saddle point are introduced in 
Sec. |TTJ In Sec. IIIII we present the results for the reso- 
nances. The calculation is performed for different exter- 
nal fields, which are varied such that lines in the two- 
dimensional space of the field strengths are traversed. 
This procedure will uncover close relations in the ener- 
gies calculated with both methods. Finally, conclusions 
are drawn in Sec. IIV1 



A. Quadratic approximation 

The transition state is located in the vicinity of the 
Stark saddle point, where the net electric force vanishes. 
For the potential caused by the external electric field and 
the nucleus, 



V t = -\+fx 



(2) 



the saddle point is located on the x-axis at x s = — l/VJ, 
and the saddle point energy has the value Vf(r s ) — 
—2^/7- The simplest approximation of the potential V) 
which contains the structure of the saddle, is a quadratic 
expansion. This expansion leads to unbound solutions 
called "Quasi-Penning" resonances by Clark et al. 0] due 
to a formal similarity with the motion in a Penning trap. 
For the readers' convenience, we recapitulate the essen- 
tial results. The expansion of the potential up to second 
order terms reads 

V f (r) = - 2y /f-y/fe+\y/f ( y 2 + z 2 ) + 0((r-r s ) 3 ) , 

(3) 

where £ = x — x s . The saddle structure of the poten- 
tial has been visualized in Fig. [TJ where the potential 
([3]) is plotted in the (£,y)-plane. The comparison with 
the exact electric potential © shows clearly that the ap- 
proximation is only valid close to the saddle point and 
can only describe states in its vicinity correctly. 

With the gauge A — (—jy, 0,0) the Hamiltonian in 
the vicinity of r s reads 

H = \ (pf + pI + pI) - ivPi + \~i 2 y 2 

-2V7+^V7V + - 2 -2C 2 ) , (4) 

and the eigenvalues of the quadratic potential yield the 
energy levels 



II. HAMILTONIAN AND RESONANCES IN 
THE VICINITY OF THE STARK SADDLE POINT 

The Hamiltonian of the hydrogen atom in crossed elec- 
tric and magnetic fields is given, in atomic units, by 



H = \{p + Af--+fx, 
2 r 



(1) 



where the electric field / is orientated along the ir-axis 
and A represents the vector potential. The magnetic field 
is supposed to be oriented along the z-axis. In this case, 
usually the symmetric gauge A = — l/2(r xje z ) with the 
magnetic flux density 7 is used, whereas, for the approxi- 
mated states discussed in this paper the application of the 
Landau gauge A = (—jy, 0, 0) is more suitable. Besides 
the energy the parity with respect to the (z = 0)-plane 
is a good quantum number, which opens the possibility 
to consider states with even and odd z-parity separately. 



En z ,ni,n 2 = ~2\/f + UJ z {n z + ^) 

+ ui(n x + -) + w 2 (n 2 + -) , (5) 

where lu z represents the frequency of the z-motion and 
is given by u z — / 3 / 4 . The separation of the coupled 
equations in the £- and y-directions with the help of an 
adequate canonical transformation [7[, leads to one real 
oscillation frequency, 

s \ ^ V2 



^ = UW-V1 +\/(7 2 -V / 7 ) +8/ 3 



and one imaginary decay rate, 



(G) 



CJ 2 



1/2 



7 2 + \/ 7 2 - V7 +8/ 



(7) 
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which describes the resonance character of the eigen- 
states. Thus, the resonance energies (real part of the 
complex eigenvalues) are completely determined by the 
two quantum numbers n z , n\ and read 

He(S) n ., ni =-2v7 + w,(n« + i)+Wi(ni + i) . (8) 

B. Exact quantum calculations 

The aim of this paper is to present evidence for the 
occurrence of quantum resonances corresponding to the 
classical transition state. For this purpose the approx- 
imate energy values ([5]) are compared with the results 
of exact quantum calculations. To perform the quantum 
calculations the Hamiltonian |T]) is rewritten in dilated 
semiparabolic coordinates 

A* = rVr + z , v = \^r - z , <p=-, (9) 
box 

where b is a convergence parameter. The transformed 
Schrodinger equation reads 

- -U 8 7 V// 2 (/i 2 + v 2 ) - 2b 6 ffiv ( f i 2 + v 2 ) cos^j^ 
= {-46 2 + A( M 2 + ^ 2 )}^ (10) 

with 

1 d d 1 d 2 

and the generalized eigenvalues A = —(1 + 2b 4 E), which 
are related to the energies E of the quantum states. We 
use a matrix representation of the full Schrodinger equa- 
tion ([10p in a complete discrete set of basis functions of 
the form 

\N^N v ,m) = \N ll ,m)®\N tl ,m) , (12) 

where \N g m) are the eigenstates of the two-dimensional 
harmonic oscillator. Resonances are uncovered with the 
complex rotation method [111, [12] • Here, the coordinates 
r of the system are replaced with complex rotated ones 
re 1 *, which leads to a complex symmetric matrix repre- 
sentation of the Hamiltonian. The complex rotation can 
be inserted in the transformed Schrodinger Eq. (|10p with 
the help of the complex extended convergence parameter 
b 

b = \b\e 1§ / 2 . (13) 

Resonances appear as complex eigenvalues E, where the 
real part of E represents the energy and the imaginary 
part is related to the width V = — 2Im(_E). 



With this method, the resonances of the hydrogen 
atom can be calculated with a high precision. It must 
be mentioned, however, that techniques for the matrix 
diagonalization used so far may fail in the energy region 
of interest. Former approaches to calculating the exact 
quantum resonances were mainly based on complex ex- 
tensions of very efficient diagonalization algorithms for 
real symmetric matrices (see, e.g., Refs. [10|, [l3() which 
require an orthogonality relation between all eigenstates. 
Of course, an orthogonality relation is not assured for 
non-Hcrmitian Hamiltonians and is, in particular, vio- 
lated at so-called exceptional points, whose existence in 
the relevant part of the spectrum has been proven re- 
cently Exceptional points are branch point sin- 
gularities of the resonances at which both the complex 
energies and the wave functions of two states become 
identical. In our calculations, the diagonalization was 
performed by applying the implicitly restarted Arnoldi 
method [l5j], which solves large scale sparse eigenvalue 
problems efficiently even for non-Hcrmitian Hamiltoni- 
ans. The matrices were built up with a basis of about 
10,000 states. 

III. RESULTS AND DISCUSSION 

To decide whether or not there is a connection be- 
tween the classical motion around the saddle point and 
features in the exact quantum spectra, we compare the 
energy eigenvalues (JSJ with the exact quantum reso- 
nances. A coincidental agreement between two energies 
of both approaches for a single parameter value is al- 
ways possible and does not answer our question. How- 
ever, a connection between the approximation (jSJ) and 
the quantum resonances is certainly present if the two 
results agree in a larger region of the parameter space. 
To check the connection we perform the calculations on 
lines (/ — fo)/l — const in the two-dimensional param- 
eter space spanned by the two field strengths, i.e., we 
introduce a parameter a and consider the field strengths 
f(a) and 7(a) as functions of a. 

An example for a comparison is presented in Fig. [2](a) 
for states with even z-parity. Here, the resonances are 
calculated for different external field strengths on the line 
defined by 

7 = 0.008 x a , (14a) 
/ = 0.0003 x a , (14b) 
< a < 1 , (14c) 

i.e., the ratio 7// of the two field strengths is always con- 
stant. The solid blue lines represent the energy eigenval- 
ues ([5]) of the quadratic approximation around the saddle 
point, and the red points are the numerically exact reso- 
nances of the Hamiltonian ([T]) . The results show that for 
the lowest quantum numbers n z , rii of the approximation 
there are very good agreements between both approaches 
on a large range of parameters a. The "Quasi-Penning" 
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FIG. 2: (Color online) (a) Comparison of the quantized en- 
ergies of the electron motion near the transition state (solid 
blue lines) with the exact quantum energies (red points) on 
the line 7 = 0.008 Xa, f = 0.0003 x a for < a < 1 as defined 
in Eqs. (|I4ap - (|I4c[l in the (7,/)-space for even z-parity. Only 
exact quantum energies with |Im(£)| < 0.0006 are shown. 
The exact resonances which trace the approximate energies 
© are shown by a larger line strength. The lowest (dashed 
green) line represents the saddle point energy. The black 
frame marks the region magnified in Fig. [3] (b) Same sit- 
uation for states with odd z-parity. 



resonances are traced by the quantum mechanically ex- 
act states when a is changed. Wherever this behavior is 
present, it is highlighted by a larger line thickness for the 
exact resonances. 

Of course, the effect does also appear for odd z-parity, 
which is shown separately in Fig.[5fb) for the same phys- 
ical parameters as in Fig.[2ja). As one can see, the agree- 
ment of the resonances is again very good in large regions 
of the parameter space. 

There are, however, parameter values a for which a 
different behavior of the energy values exists. For low 
field strengths the Coulomb potential becomes too domi- 
nant, and the expansion of the potential in the vicinity of 
the Stark saddle point no longer leads to reasonable solu- 
tions, as can be observed in Fig. [2] Furthermore, avoided 
level crossings of the numerically exact resonances have 
no counterpart in the approximative states. Evidently, 
couplings of states, which are the origin of avoided level 
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FIG. 3: (Color online) Same situation as depicted in Fig. [2] 
The region around the avoided level crossings near the state 
with n z = 2, m — is magnified. 




FIG. 4: (Color online) Same comparison of the quantized 
energies of the electron motion near the transition state (solid 
blue lines) with exact quantum energies (red points) as in Fig. 
[2] but on the line 7 = 0.01 x a', f = 0.0001 + 0.0002 x a' for 
< a' < 1 defined in Eqs. (|15ap - (|I 5c|l . Only exact quantum 
energies with |Im(£')| < 0.0002 are shown. The saddle point 
energy is marked by the lowest (dashed green) line. 



crossings, are not included in the simple model. Thus, 
the deviations in these regions are not surprising. Fig. 
[3] shows an example of a region with avoided level cross- 
ings of quantum mechanically exact resonances. It can 
be clearly seen that far from the crossing the agreement 
between the numerically exact energies and the energy 
value corresponding to n z = 2, n\ = in Eq. ([8]) is very 
good, whereas at the avoided crossings themselves the 
approximative solutions cannot reproduce the results of 
the full Hamiltonian |T]). 

The calculations shown in Fig. 3] are performed on the 



5 



TABLE I 


: Quantum 


numbers 


n z and ni 


of the "Quasi- 


Penning" 


resonances 


closely related with 


exact quantum 


states. 










n z 




ni 


n z 


ni 










4 





1 










1 


2 










2 


3 












line 

7 = 0.01 x a' , (15a) 

/ = 0.0001 + 0.0002 x a' , (15b) 

< a' < 1 , (15c) 

i.e., the ratio of the field strengths is no longer constant. 
For a' — there is still an electric field present but no 
magnetic field. The important property of this line is the 
existence of the saddle point for all values of a' drawn in 
Fig. 21 Again, there are exact quantum resonances which 
behave like the approximate solutions around the Stark 
saddle point and lead to a good agreement between both 
methods on a wide range of parameter values. Also near 
a' = 0, where a pure Stark effect is present, the classical 
electron motion near the transition state describes the 
quantum resonances very well. 

The calculations performed here reveal that a close 
connection between the two approaches exists. It ap- 
pears for the seven pairs of quantum numbers of the 
"Quasi Penning" resonances listed in Table|H This rather 
high number of concordant resonances shows that, in- 
deed, there is a signature of the classical transition state 
in the exact quantum spectrum. 

IV. CONCLUSION 

In summary we have demonstrated that there is clear 
evidence for the correspondence of classical motion of 



electrons confined to the vicinity of the Stark saddle point 
and resonances in the exact quantum spectra. The agree- 
ment of the energies in the classical transition state ap- 
proach and the quantum calculations is very good. This 
appears surprising because the simplest expansion of the 
potential, only up to second order, was used, i.e., the ap- 
proximation is only valid in a small region very close to 
the saddle point. This clearly can be understood as a sign 
for an eminent impact of the transition state on the exact 
quantum spectrum, and demonstrates the importance of 
the transition state for the ionization mechanism. 

It must be noted, however, that the excellent agree- 
ment between both methods exists only for the real parts 
of the complex energy eigenvalues. The simple approxi- 
mation used in the expansion of the potential is not ca- 
pable of reproducing the correct decay rates or imaginary 
parts. 

As an outlook, it seems worthwhile to extend the calcu- 
lations to the normal form procedure developed by Uzer 
et al. @ and to check whether higher orders can lead to 
an even better description of the quantum resonances by 
classical electronic motion. As can be expected, only the 
lowest "Quasi-Penning" resonances turn out to have a 
counterpart in the exact quantum resonances. It will be 
interesting to investigate which are the highest quantized 
levels of the transition state theory with a correspondence 
in the quantum spectrum. 
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